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Abstract 

In quantum computation, it is of paramount importance to locate 
the parameter space where maximal coherence can be preserved in the 
qubit system. Considerable insight in environment-induced decoher- 
ence has been gained in the last decade from detailed studies using the 
quantum Brownian motion (QBM) models. A number of respectable 
authors have applied these insights derived from QBM models to two 
level systems interacting with a field. Their conclusions based on 
this particular type of qubit coupling to the environment had led to 
the general belief that 2LS are easily decohered. In a recent paper 
, we debunk such a myth and caution indiscriminate application of 
the QBM model of decoherence to arbitrary two level systems. We 
point out that at least for a two-level atom (2LA)- electromagnetic 
field (EMF) system alone, as used in the atom cavity prototypes of 
quantum computers, the decoherence time is rather long, compara- 
ble to the relaxation time. In the standard Hamiltonian of the 2LA, 
the dominant interaction is the a± type of coupling between the two 
levels (what constitutes the qubit) and the field, not the (T^ type as- 
sumed in most previous discussions of qubit decoherence, which shows 
the QBM behavior. Depending on the coupling the field can act as 
a resonator (in an atom cavity) or as a bath (in QBM) and produce 
very different decoherent behavior. Our conclusion is based on a new 
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exact master equation we derived at zero temperature which general- 
izes the text-book ones restricted by the Born-Markov approximation. 
Indeed many cavity experiments testify to the correctness of these re- 
sults, and that the 2LA-EMF system maintaining its coherence in 
sufficiently long duration is the reason why experimentalists can ma- 
nipulate them to show interesting quantum coherence effects. 



1 Introduction 



A two-level system (2LS) interacting with an electromagnetic field (EMF) has 
proven to be a very useful model for a wide range of problems from atomic- 
optical [0, ^, |, H, 0, H and condensed matter |10| processes to quantum 
computation For the latter application stringent limits in maintaining 
the coherence of the 2LS (called qubits) are required. This prompted us 
to revisit the theoretical structure of the 2LS-EMF system, paying special 
attention to its relaxation and decoherence properties [|12|, [l^, |l^ |n 

m 
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Environment-induced decoherence has been studied extensively in re- 
cent years primarily based on models of quantum Brownian motion (QBM) 

l8i for the interaction of a simple harmonic oscil- 
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23, |23, g|, gl, g|, 

lator (Brownian particle) with a harmonic oscillator bath (HOB) at a finite 
temperature, leading to a reasonably good understanding of its characteristic 
features. Decoherence of a 2LS in an EMF has been studied by a number 
of authors, notably |jl2|, |13|, [l^ , and their dissipative and decoherent behav- 
ior are reported to be similar to that of a QBM in a harmonic oscillator 
bath. They described the progression in three stages - quiescent, vacuum 
fluctuation-dominated and thermal fluctuation-dominated, separated by the 
cutoff frequency and the thermal de Broglie frequency (wavelength), which 

l2i. A recent com- 



are characteristic QBM features fS^, ^ 



prehensive review of environment-induced decoherence can be found in [53 

Our findings, based on the standard 2LA-EMF model 0, are in stark 
disagreement from that reported in the literature. From solutions of the 
exact master equation for the reduced density matrix obtained recently by 
Charis Anastopoulos and the author |]l| capable of treating non-Markovian 
dynamics at zero temperature, the decoherence time of the 2LA is found to 
be close to the relaxation time. In many circumstances these times are rather 



2 



long - indeed, that is why observations of many optical coherence phenom- 
ena are possible. This result is, in first appearance, rather counter-intuitive, 
and different from all previous findings. However, as will be explained here, 
the 'intuition' one forms about the nature of dissipation and decoherence 
are hitherto largely based on the QBM model. It is quite general, yet by 
no means universal. This convenient extrapolation could have influenced the 
choice of model in previous authors' investigation and the readers' impression 
of decoherence for a 2LS. Our findings show that such a commonly invoked 
intuition for QBM in a HOB fails to apply to that of a two-level atom (2LA) 
interacting with an electromagnetic field (EMF). Our conclusion shows the 
need for extra caution in accepting categorical depiction of decoherence: One 
needs to identify the correct type of coupling of the 2LS with its environment 
and consider different decoherent mechanisms at work in a realistic experi- 
ment. 



Decoherence in QBM 



The folklore that decoherence of a Brownian oscillator proceeds in a very 
short time (typically 10~^° of the dissipation time) is really based on high 
temperature ohmic bath conditions. This is widely known because it is the 
only technically simple case studied in detail |Q It is by no means uni- 
versal. Intuitively, the bath needs to have many degrees of freedom, prefer- 
ably acting independently of each other so incooperatively that the phase 
information in the system will be dispersed to the largest extent amongst 
the many bath degrees of freedom and affords little chance or takes inordi- 
nately long time to be revived or reconstituted (recoherence, see, e.g., p^ ). 
Long decoherence times can appear in cases of low temperature, supraohmic 
bath, as was first pointed out by Hu, Paz and Zhang |^in analyzing deco- 
herence behavior of QBM based on their exact non-Markovian master (HPZ) 
equation for general environments. (For details based on numerical and an- 
alytical solutions of the HPZ equation, see p9|, ^). The almost opposite 
picture (of very long decoherence time) is exemplified by two coupled sub- 
systems where no coarse-graining is introduced, such as the 'dressed atom' 
description made possible because of the coherence established between a 
2LA and a single resonant mode in the cavity. (Collapse and revival of the 
Rabi nutation are distinct features of quantum coherence 0.) This shows 
that both the coupling and the nature of the environment are important 



3 



factors affecting the decoherence of the system. 

We also want to point out that equating a 2LS with QBM could be as 
mistaken as ignoring the basic difference between dissipative and mixing sys- 
tems. To the extent that QBM exemplifies the former, spin echo phenomena 
(e.g, Chap. 3 [§])is an example of the latter. Plasma waves showing Lan- 
dau damping in Vlasov dynamics arising from a mean field approximation 



is another example [pSj |36[| of the latter - the 'damping' is really a mis- 
nomer. Just as in the spin echo phenomena, the basic physics in this case is 
not dissipation in the Boltzmann sense, but statistical mixing We will 



see that the statistical mechanical properties of a 2LA-EMF system is closer 
to the latter than the QBM, and reflects in their different decoherence nature. 

Coherence in the 2LS 

For the 2LA-EMF system, one clear distinction between an EM field and 
a system of harmonic oscillators acting as bath is that the field (coupled to a 
detector) has an intrinsic spectral density function, which cannot be chosen 



arbitrarily. For example, it has been shown |2g] that a conformal scalar field 
in two dimensions coupled to a monopole detector has an Ohmic character 
while in four dimensions it is supraohmic . Barone and Caldeira |^ showed 



that the spectral density function for EM fields with momentum coupling 
to an oscillator detector is supraohmic. These density functions would show 
very different decoherent behavior from the high temperature Ohmic HOB 
case which molds the common impression of decoherence. 

But the most important distinction from QBM is that the 2LA couples 
with the EMF in the discrete number basis for the field, unlike the continuous 
amplitude basis in the QBM. This fact (which is true in the rotating wave and 
dipole approximation) implies that the 2LS plus EMF system is a resonant 
one. Hence even though the EM field has just as many (in fact, a large 
number of) modes as the HOB, only a very small fraction of them in a 
narrow range of the resonance frequency are efficiently coupled to the atom. 
This is the root cause for the very different qualitative behavior between the 
QBM and the 2LS as far as decoherence is concerned. 

We see that in a 2LA, for purely radiative decay the decay time Ti of the 
inversion is half the decay time T2 of the polarization. There is no large order 
of magnitude differences between dissipation and decoherence time. In fact 
it is perhaps inappropriate to talk about dissipation for a 2LA-EMF system 
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because the conditions for a bath to actuate such a process is lacking. The 
transition from excited to ground state is closer in nature to relaxation (with 
relaxation time constant F) than dissipation. In a cavity where excitation of 
the atom from the field (absorption)balances with emission, it is more appro- 
priate to refer to the resonant state of the atom-field as a coherent system. 
In these scenarios the distinction between QBM and 2LA cannot be clearer. 



Difference between QBM and 2LS 



So what led earlier authors to make the claim that 2LS decoheres easily? 
We think the confusion arises when the picture of QBM dissipation and 
decoherence is grafted on the 2LA-EMF system indiscriminately. If the field 
which acts as the environment is a phonon field (from ion vibrations, see, 
e.g., |TB[), if there is atomic collisions in a cavity H] or if the cavity walls are 



imperfect, some of these decoherent agents could follow the QBM pattern as 
reported by many authors. Such sources can be important for some setups. 

Quantitatively, the model for the 2LS used by most authors for the dis- 
cussion of decoherence inspired by QBM type of behavior has the atom in 
a dz state (the diagonal Pauli matrix) coupled to the field mode operators 
¥ ,h. This type of coupling term (call it type for convenience) commutes 
with the Hamiltonian of the system, and admits a diagonalization in the 
eigenbasis of the Hamiltonian. The field is coupled to the atom as a whole, 
not to the two level system (qubit), which is at the heart of quantum com- 
putation. By contrast the standard model for 2LA-EMF which we studied 
has a a± coupling to the field modes which controls the two- level activity of 
the atom. This coupling considered in the standard model is indispensable, 
i.e., it cannot be removed from the two- level atom as it defines it and will 
be present in any realistic situation. Where then could the QBM type of 
interaction enter in the 2LA? 

If the EM field is the only environment present, we can still ask if a QBM 
type of coupling term with the EM field would appear, and if yes, how strong 
would its effect be? A useful way to compare the relative importance of these 
two types of coupling would be to seek out the source where they stem from. 

Recall that the standard model is derived under the dipole and rotating 
wave approximations. In the next section we will show that the type of 
coupling appears only in the next order expansion after the dipole approx- 
imation. Since these are good approximations for a large class of atomic 
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states when the atom is nonrelativistic, the contribution from the QBM type 
of couphng used in [|l^, |13|, |l^, [l^ should be neghgible in a 2LA-EMF sys- 
tem and its ensuing decoherent effect insignificant. In this sense the EM field 
does not in leading order of approximation act like a bath in the QBM way, 
and coherence in a 2LA-EMF system is quite well preserved (excepting other 
processes, e.g. 



T|,|T1). 



Decoherent behavior of the 2LA-EMF system we are reporting on here is 
based on an analysis of the full (non-Markovian) dynamics in the cases of a 
free quantum field and a cavity field at zero temperature from solutions of an 
exact master equations derived by Anastopoulos and Hu They used the 
influence functional method 0] to take into account the full backreaction of 
the field on the atom, while adopting Grassmannian variables for the 2LA 
and the coherent state representation for the EMF. 



2 The Model 

Our model for atom-field interaction is the standard one (see Appendix A 
of for details) 0, ^ |[ [|. The total Hamiltonian for a (stationary) atom 
interacting with a quantum electromagnetic field (EMF) under the dipole, 
rotating wave (RW) and two-level (2L) approximation is given by 

H = hiVoSz + nY^ [iv^bib^ + {g^Sj^ + g^S Ji)] (2.1) 

k 

where fej^, 6k are the creation and annihilation operators for the kth normal 
mode with frequency cuk of the electromagnetic field (thus for the field vacuum 
^k|0) = 0, [&k5 ^k'] — '^k,k'5 for all k.), and ujq = U21 is the frequency between 
the two levels. Here 

1 ^ 1 
Sz = -^z, S± = cr± = -{a^ ± iay) (2.2) 

where (Tx,y,z are the standard 2x2 Pauli matrices with = diag{l, —1), etc. 
The coupling constant g^ = c?2ik/k(X) where 



^Our Hamiltonian is given in the so-called minimal coupling (MC) as different from the 
multipolar coupling (MP) 0, which may be more relevant to atoms in a cavity because 
the explicit Coulomb interaction between the atom and its image charge is removed. 
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and dij = e J (f)iX.(f)jdr'x is the dipole matrix element between the eigenfunc- 
tions (pi of the electron-field system, eko- is the unit polarization vector ( 
a = 1,2 are the two polarizations), and /k(x) is the spatial mode functions of 
the vector potential of the electromagnetic field (in free space, /k(x) = e~^^'^, 
V is the volume of space.). Under the dipole approximation is evaluated 
at the position of the atom X. Since dj^ = d^j, dij^ = dji-^, we will choose a 
mode function representation such that (/k is real. When only one mode in 
the EM field is considered, this is the Jaynes-Cummings model. 

To see how this could possibly be related to the type of coupling with 
Hamiltonian (used by e.g., ||T3|, [I^ for the study of decoherence in 2LS) 

H = huoS, + hJ2 [^k^k^k + ha, (gkbk + ^k&l)] (2-4) 

k 

we examine the next term after the dipole approximation. This has a con- 
tribution to gijk even when i = j. This is equal to 

gnk = Ckk ■ qi (2.5) 

where 

^i = Y. 'Pi^^iP ■ ek^)<Pidx^ (2.6) 
and Cfe is a constant given by 

Ck = --(2;i^keo^)~'/' (2.7) 



m 

It generates an additional coupling term 

^zigikbk + gikbi) + l(5'2k&k + 92kbi) (2.8) 

k 

where 

gik = fi'llk - fi'22k, 5'2k = gUk + 5'22k (2.9) 

This gives the lowest order a, type of coupling in a 2LA -EMF system. The 
ratio of the coupling g^^ of the a, type in (|2.4|) to the dipole coupling gi^ in 
:TDis 

I / I | k(qi-q2) | ^ ^k|qi-q2| 

W^k = 1 < 1 (2.10) 

muj]f^di2 ma;k"i2 

Thus the a, type of coupling generated from the 2LA- EMF interaction will 

be significant only for very high frequencies cuk of the EM field, a point 

intuitively clear from the meaning of the dipole approximation. 
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3 The Master Equation 



We use the open system approach and the influence functional formalism to 
derive a master equation for the 2LA with backreaction of the field treated 
self-consistently. This involves deriving the influence functional and the evo- 
lution operator for the reduced density matrix of the 2LA. For convenience 
we use Grassmannian variables for treating fermions, and the coherent state 
representation for the field. The coherent state of the combined atom-field 
system is 

\{z},v) = \{z})y<\v) (3.1) 

where \z), z a complex number, denotes the EM field coherent states and 
\r]), f] a Grassmannian or anticommuting number, denotes the electron co- 
herent state. We assume initially that the density matrix of the total sys- 
tem-|-environment is factorizable /3(0) = Pe(0) ® Pb{0). Only at that time 
would z and t] be pure complex and Grassmannian numbers respectively. As 
the system evolves, both t] and z contain Grassmann and c-number parts. 
The mixing of even and odd parts (note g^^ is odd) comes about as the initially 
factorized atom state becomes " dressed" . 

We skip over the details of the derivation which can be found in the orig- 
inal paper but simply present the master equation, here in operator form 
(at zero temperature with the field in a vacuum state) as follows: Writing 

^ = r(t) + zm (3.2) 
u{t) 

the master equation reads 

= -t[Hit),p] + r(t){5+5_,p} - 2Tit)S.pS+ (3.3) 

where 



H{t) = Q{t)S+S^ (3.4) 

The first term corresponds to the unitary Hamiltonian evolution, only now 
the effect of the environment has induced a time dependent shift in the value 
of the frequency, the second term is time dependent dissipation and the third 
corresponds to noise. 
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The effect of the field is contained in the u{s),u{s) functions, which are 
obtained as solutions of the linear integro-differential equations 

■f] + iuri + / ds'fi{s — s')r]{s') = (3.5) 
Jo 

fl'-iujf]'+ r ds'fj.*{s- s')f]'{s') = (3.6) 



under the condition 

u{0) = u{t) = 1 (3.7) 
Here, the kernel fi{s) is given by 

/^(s) = E^?ke-^^''^ (3.8) 

k 

The equations for u{t) functions can be solved with the use of the Laplace 
transform and the convolution theorem. They are given by 

uis)=C-^ : — - =— -/ ; — - 3.9 

\Z + lUJ + IX{Z) ) ZTTt Jc-ioo Z + tUJ + fi{zj 

where fl{z) is the Laplace transform of the kernel fi{s) and c is a real constant 
larger than the real part of the poles of the integrand. 

3.1 Spontaneous emission 

To show how the standard results are regained, and to understand the mean- 
ing of the new function in the master equation, let us consider the physical 
process of spontaneous emission. Start with a generic initial density matrix 

its corresponding Q-symbol is 

p(f/, T]) = X + y*ri + yfj + {1 — x)f]r] (3-11) 

If we evolve it with the density matrix propagator (derived in []I|) we obtain 
for the state at time t 

Ptiv, v) = ^ ~ uu{l — x) + {uy*ri) + {ufjy) + (^^(1 — x)) fji] (3.12) 
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corresponding to 



/ uu{l-x) uy \ 

uy* l-uu{l-x) ) ^"^-^^^ 

Considering the case x = ?/ = we get for the probabihty of spontaneous 
emission 

P(l ^ 0,t) = 1 -MM (3.14) 

The rate of decoherence in the energy eigenstates is governed by the absolute 
value of the function u (the off- diagonal terms) while u itself determines the 
rate of energy flow from the atom to the environment. Hence for our partic- 
ular choice of initial state (vacuum) we find that decoherence and relaxation 
time are essentially identical . This equation is useful for studying decoher- 
ence of a qubit in a QED cavity. 



3.2 Field modes in Free Space and Cavity 

Our master equation ( p.3|) depends solely on the function u{t), which in 
turn is determined by the kernel /i(s). In [|I| we have given some analytic 
expressions for this function in various cases, including a single mode, an 
infinite number of modes in the field, and a cavity consisting of two parallel 
plates at distance L. The field satisfies Dirichlet boundary conditions on the 
surface of the plates. In all cases our results obtained from the generalized 
master equation bear closer resemblance qualitatively to the 2LA behavior 
than that of quantum Brownian motion. Details can be found in |jl|. We 
now summarize the major findings obtained there. 



4 Discussion 

The physics of a 2LA-EMF system at zero temperature is characterized by 
a number of time constants: 

1) The inverse natural frequency uj^^ 

2) The inverse coupling constant g^^ = ^/uju./ A 

3) The relaxation time constant F"^ 

4) The cavity size L (divided by c) 

First consider a zero temperature field in free space, thus ignoring factor 
4). Start with only one mode in the field in resonance with the atom (whence 
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we omit the subscript in (yf^), then the system undergoes Rabi nutation with 
frequency Q ~ ^fi/rT+T, where n is the photon number in the field. The 
collapse time (assuming a large mean photon number n) is g~^, and revival 
time is 2-n^/fi/g. §. At om excitation becomes significant in a time much 
greater than uJq^ but shorter than g~^. (This is the condition for a first order 
perturbation theory to give reasonable results.) For a large number of modes, 
spontaneous emission occurs at the relaxation time scale = 1 9 » t^o"^ 
which we found to be the same as the decoherence time - the time for the off 
diagonal elements of the reduced density matrix to decay. When the mean 
number of photons in the field is large {n » 1), they become comparable 
to the collapse time. This is a measure of the coherence in the atom-field 
system, and is controlled mainly by their coupling and the photon number in 
the field. We see that with the resonance condition, the nature of decoherence 
in 2LS is very different from the QBM situation, where phase information in 
the Brownian particle is efficiently dispersed in the many modes in the bath 
coupled almost equally to the system. As we remarked in the Introduction, 
the identification of the phase information and energy fiow from the 2LS to its 
environment is similar to the spin echo phenomena (Landau 'damping') which 
is based on statistical mixing rather than dissipation. The mathematical 
distinction lies between considering the system coupled to the continuous 
amplitude basis (QBM) of the environment or to the discrete number basis 
(2LA). The QBM case essentially produces noise that drives the system in 
a way insensitive to its own intrinsic dynamics. In our model, the coupling 
respects the internal dynamical structure of the 2LS and allows it to keep its 
coherence. 

To see how the distribution of modes in a field changes the picture, the 
cavity field calculation is useful. There, as the plots in show, the relax- 
ation constant develops peaks and minima. The resonance effect is enhanced 
by a cavity size commensurate with the natural frequency of the 2LA and dis- 
sipation weakens. Narrow band resonance fiuorescence as well as inhibition 
of spontaneous decay by frequent measurements - the Quantum Zeno effect 
- are interesting phenomena which our equations can provide finer details. 

Non-Markovian processes involve memory effects (nonlocal in time). When 
the reaction time of the bath is comparable to or faster than the natural 
time scale of the system {ujq), one also expects to see non-Markovian be- 
havior. For the QBM problem, the case of high temperature Ohmic bath 
is almost the only condition that would yield a Markovian dynamics. For 
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other types of spectral density (supraohmic) or for a low temperature bath, 
the dynamics of the system is generally non-Markovian |2^. By contrast, 



the 2LA is quite different: At zero temperature there is only one time scale 
= X~'^uj~^ » uj~^ that determines both decoherence and relaxation. 
There is no memory effect and hence the process is Markovian. We expect 
that at finite temperature the dynamics of the 2LA will be nonMarkovian 
This is because there are more ways for the atom and the field to 
get entangled, and the memory effects of their interaction would presumably 
persist. 

To conclude, we show that not all two level systems (2LS) decohere hke 
quantum Brownian particles (QBM). Specifically, in the case of a two- level 
atom (2LA) interacting with an electromagnetic field (EMF) alone through 
the standard a± coupling as in a quantum optics / cavity QED quantum com- 
puter prototype, decoherence time of the 2LA is close to the relaxation time. 
This behavior is completely different from the QBM. The crucial difference 
lies in the type of coupling between the 2LS and the environment: The 2LA 
interacts resonately with selected modes in the electromagnetic field while 
the QBM is coupled to all modes in the harmonic oscillator bath. Only for 
those mechanisms which can be described by a type of coupling of the 2LS 
with the environment (which is subdominant in a 2LA-EMF system) will one 
expect a QBM-like decoherent behavior. Therefore one has to be careful in 
specifying what type of 2LS constitutes the qubit and and how it couples to 
the field, since the field can act as a resonator (2LA) or as a bath (QBM) 
producing very different decoherent behavior. In a realistic experiment one 
needs to know which processes are describable by what type of coupling and 
the weight of their relative contributions in varying environmental conditions 
before drawing credible conclusions about the overall nature and degree of 
decoherence in that system. 
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